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A complete thermodynamic treatment of the Casimir effect is presented. Explicit ex- 
pressions for the free and the internal energy, the entropy and the pressure are discussed. 
t^j- ' As an example we consider the Casimir effect with different temperatures between the 

plates (X) resp. outside of them (X"). For X" < X the pressure of heat radiation can 
eventually compensate the Casimir force and the total pressure can vanish. We consider 
both an isothermal and an adiabatic treatment of the interior region. The equilibrium 
point (vanishing pressure) turns out instable in the isothermal case. In the adiabatic 



situation we have both an instable and a stable equilibrium point, if X'/X is sufficiently 
L^ ■ small. Quantitative aspects are briefly discussed. 

^ PACS: 05.90.+m, 11.10. Wx, 61.16.Ch 
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1. Introduction 

The Casimir effect [|TJ is one of the fundamental effects of Quantum Field Theory. It 
tests the relevance of the zero point energy: two parallel, large conducting plates in 
a distance a change the vacuum energy of Quantum Electrodynamics in such a way, 
that a net attractive force between the plates results. Qualitatively, this situation 
does not change, if the system is at finite temperature. The quantitative corrections 
for finite temperature [§], finite conductivity [|3| and the necessary changes for related 
problems (e.g. for different geometries) [[|] are known, and the Casimir force has been 
experimentally established || || [/J. The thermodynamics for the original setup has 
been investigated previously P, || ||. Here we shall give a somewhat more complete 
treatment, in which we consider the region between the plates and outside of them 
separately. This allows for a situation, in which the two regions have different temper- 
atures (outside T", inside T). If we take T' < T, the external pressure is reduced in 
comparison with the standard situation (T" = T). Therefore we expect the existence 
of a certain distance a , at which the Casimir attraction is compensated by the net 
radiation pressure. We shall investigate this (mechanical) equilibrium point and its 
stability both for an isothermal and an adiabatic treatment of the interior region. 

2. Thermodynamic functions 

In this section we shall collect and discuss briefly formulae for thermodynamic functions 
for the original Casimir setup. A major part of the results has been found before J8], [| by 
various methods. In order to make this paper self-contained, we shall indicate briefly 
a standard procedure for the derivation and list the results. In addition, we shall 
give numerical details on the functions involved. We consider two parallel, perfectly 
conducting square plates (side L, distance a, L > a), embedded in a large cube (side 
L) with one of the plates at face and periodic boundary conditions. We shall consider 
contributions from the volume L 2 a between the plates (suffix int) resp. L 2 (L — a) 
outside of them (suffix ext) separately. In terms of the partition function Z the free 
energy of a photon gas at temperature T reads 

F = ~lnZ = Y,(Go(k) + G p (k)) (1) 

with 

he l 

G {k) = —k, G p {k) = - In (1 - exp (-/3hck)) . (2) 



Here (3 = 1/ksT with Boltzmann's constant ks and tick = hcvk 2 is the energy of 
a photon with wave number k. The sum extends on all accessible wave numbers and 
polarisation states. In order to keep track of ultraviolet divergencies, we shall regularize 
the contribution at temperature zero replacing k — > A; exp (— Xk) in Go and considering 
A — » at the end. We assume L fixed, but so large, that sums on photon wave numbers 
can be replaced by corresponding integrals. Then the contributions to the free energy 
per unit area <fi = F/L 2 read. 
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The integrals can be done analytically. In the results to be listed below we have 
expanded the result for Go in powers of A and omitted all terms, which vanish for A — > 0. 
For G/3 it turned out convenient to use the polylogarithm C n {z) = J2m=i z m /m n , which 
allows even for a calculation of the indefinite integrals and has simple properties (cf. 
Appendix). Further thermodynamic functions can be calculated from by standard 
formulae. We shall give results for the entropy 
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the pressure 



and the internal energy 
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The results are listed in Table 1 in terms of a finite function g(v) to be specified below. 
The variable v is the dimensionless quantity (cf. p, |9|, [IT 

a aks 



iihcj3 nhc 



T. 



(8) 



Table 1 
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The ultraviolet behavior can be discussed without any details on g. Both contributions 
to a are finite. The divergent contributions to the pressure cancel in the sum -P ex t + -Pmt 
and in the difference P ex t{P') ~ -fext(/5), which we shall use in the next section. The 
total free energy ext + int contains a divergent constant ~ L/X A , as well as the total 
internal energy. 

Another formula, which can be read off Table 1 without knowledge of g is the relation 
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which is connected with the fact, that physically relevant contributions scale as h(v)/a? 
with some function b. 



The function g(v) contains an infinite sum, which cannot be performed analytically. 
Various forms can be found in the literature || ||. We shall prefer a representation, 
which allows for a fast and accurate numerical computation. For this purpose we use 
the functions (cf. Appendix) 
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h(x) = xk'(x) 
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For x > the sum converges rapidly, k is positive and h is negative. Both functions 
vanish exponentially for large argument. Two representations of g in terms of k read 
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with C(3) = £3(1) = 1.2020569 • • •. The two forms are related by Poisson's sum formula 
[fP2| . Taken together, they contain information on the behavior of thermodynamic 
functions under temperature inversion (T — > 1/T) || |l'l] (the quantities used in || are 

£ = ™, /(0 = vrV/45 + 7r 2 (7(t;)). 

In order to discuss the behavior of the thermodynamic functions quantitatively, we 
write 
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<P = <p L (T) + —-f(v), a = o L (T) + -s(v), P 
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p(v), e = e L (T) + 
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(13) 
where (4>l, &l, &l) refer to the extensive contributions ~ L from the external cube. The 
remainders are connected by simple relations. As a consequence of (9) we have 



3f(v ) + v s(v ) — p(v) 
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(14) 



Equ. (7) amounts to 

e{v)= f{v)+vs{v). (15) 

Therefore only two of the four functions (f,s,p,e) are linearly independent. Explicit 
forms in terms of g can be read off Table 1. Expressions in terms of k(x), h(x) result by 
insertion of g. There are two equivalent forms (A,B) for every function corresponding 
to (11) resp. (12), which are listed in Table 2. 
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The leading terms for small values of v (i.e. at low temperatures) are found from forms 
(A) neglecting (k, h). This approximation can be used for v < 0.09 with errors less than 



1%. The asymptotic behaviour at large values of v (i.e. at high temperatures) is ob- 
tained from forms (B) neglecting (k, h). The corresponding asymptotic approximations 
for (/, s,p) can be used in a rather large domain. For v = 0.25 the errors are <1% and 
decrease rapidly with growing v. Form (A) shows, that the entropy vanishes for T = 0, 
i.e. Nernst's law is fulfilled. On the other hand T tends to a finite value T(a = 0) for 
a — > 0, which is determined by the entropy alone. Form (B) shows, that the entropy 
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It has to 



becomes constant in the high-temperature limit (Kirchhoff's theorem 
be noted, that there is only a rather narrow domain in v, in which the infinite sums 
(k, h) play a role. In this region the low-temperature behavior of the Casimir force is 
gradually changed by temperature effects. This domain is realistic || [7|: for T = 18°C 
we have v ~ 0.04a(/im). Using the expression (10), one can find numerical values in 
this domain rapidly and with high precision. Results are shown in the plots given in 
Fig.l (in which the broken lines correspond to the asymptotic approximations). 
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Fig.l Thermodynamic functions. 



3. Different Temperatures and Equilibrium 



As an application we shall now investigate a situation, in which the temperature T 
between the plates may differ from the temperature X" of the exterior region. We 
shall keep T" resp. v' = aksT' /nhc fixed. For the interior region we consider two 
possibilities. In the first case also T resp. v is kept fixed (isothermal situation). In 
the second case we treat the internal region as a closed system, so that the entropy 



<Ji n t remains constant (adiabatic situation) and T resp. v vary appropriately with a. 
In both cases we shall focus attention on the question, whether there is a distance do, 
at which the total pressure vanishes. If one of the two plates is movable, the distance 
a = a would constitute mechanical equilibrium between the Casimir attraction and 
the thermal radiation pressure (which repels the plates). We shall also consider the 
question, whether this equilibrium is stable. 



For an appropriate treatment we consider the quantity 
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in terms of p from section 2. We shall always use the form (B) from Table 2. The 
equilibrium condition 



P(a, v,v') 
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amounts to 
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Solutions are only possible, if v' < v(a), since p is negative. A Taylor expansion of P 

near a = a shows, that dP/da at equilibrium must be negative for a stable situation. 

Differentiating (16) and using (19), one observes, that the stability condition amounts 

to 
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We shall always seek solutions considering a scaled form of (19) 
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with some scale vq. 



In the isothermal case v is proportional to a. If we take vq = v, the ratio k = T'/T 
does not depend on a. From the plot for p(v) given in Fig. 1 we can see, that R(v) is 
a monotonous function with positive tangent. The function is positive above a certain 
value of v. For any given value of < K < 1 we have therefore one solution, which is 
unstable. Some numerical results are listed in Table 3. The third row gives the distances 
for T = 18°C in /im. The last row is obtained using the asymptotic approximation 
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-i>C(3)/47r , which turns out quite accurate. 



Table 3: isothermal equilibrium 


K 





0.2 


0.4 


0.6 


0.8 


0.95 


v(a) 


0.2419 


0.2421 


0.24340 


0.2532 


0.2879 


0.4233 


a(18°C)//im 


5.981 


5.984 


6.032 


6.260 


7.117 


10.46 
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In the adiabatic case we have to evaluate (21) along an isentrope, which corresponds 
to a curve &i n t( v ) — const. A convenient parametrisation is obtained using the constant 
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instead of <Tj ni . An equation for the isentrope is obtained inserting the explicit ex- 
pression for (Ti n t(v) = —7rg'(v)/a 2 from section 2. In terms of t and s the equation 

reads 
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The leading term on the r.h.s. for small v is v 2 (see Table 2, form (A)). Therefore the 
isentrope becomes an isotherme for small v and t can be identified as 
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Fig. 2 Isentrope. 



In order to obtain an explicit form for the isentrope we have to solve (23) for v = v(at). 
This can be done numerically with high precision. We have used the form (B) from 
Table 2 for s(v) to determine v(at). The result is plotted in the left diagram of Fig. 2. 
An approximate form v as (at) is obtained neglecting k + h in s(v). The result is 
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v can be approximated by v as for at > 0.25 with an error < 0.1%. With increasing 
at the error decreases rapidly. The right diagram of Fig. 2 (in which the broken line 
corresponds to v as ) illustrates this fact. 

If we take vq = at in (21), the ratio k(0) = T'/T(a = 0) does not depend on a and can 
therefore be used as an input for (21). The function R(v) shows, however, a different 
behavior than in the isothermal case. It vanishes for at = 0.2763, assumes a maximum 
at at = atM = 0.4391 and falls off slowly for larger argument (see Fig. 3). 
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Fig. 3 Adiabatic equilibrium function R. 

Therefore we obtain no solution, if ac(0) is larger than % = (R(atM)) 1 ^ 4 = 0.68542. For 
k < km there are two solutions: one of them (corresponding to the lower value of at) 
is unstable, the other one is stable. The distance between the two solutions increases 
with decreasing /t(0). Some numerical results are listed in Table 4. The entries for 
a(18°C) refer to T(0) = 18°C. 



Table 4: adiabatic equilibrium 




«(0) 


0.2 


0.4 


0.6 


0.65 


unstab. 


at 


0.2767 


0.2819 


0.3148 


0.3441 




a(18°C)//im 


6.839 


6.969 


7.783 


8.506 




«(a) 


0.2285 


0.4592 


0.7093 


0.7876 


stab. 


at 


26.03 


3.235 


0.8917 


0.6503 




a(18°C)///m 


643.4 


79.98 


22.05 


16.08 




n(a) 


0.99998 


0.9984 


0.9778 


0.9565 



It is observed, that n(a) is close to 1 for the stable equilibrium points. This is due to 
the slow falloff of R(v). Using the asymptotic approximation p — > — t> as (a£)C(3)/47r 2 , 
the data in Table 4 are reproduced with rather small errors. 

The existence of stable equilibrium points implies, that a movable plate may perform 
oscillations about these points. In order to approach the problem of an observation, 
we shall discuss briefly some quantitative details. For these, we have used the formula 
for the Casimir pressure 

n 2 hc - N 

1.3001 • 10~ 7 — - (26) 
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Here a is measured in microns (/mi). We assume the movable plate situated at the 
distance a\ from a stable equilibrium point a, where a\ja is small enough, so that 
a linear approximation is possible. For the data in Table 4 (stable points, 18°C) we 
obtain a restoring force/area of 



^(0.39- 28)^, 
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The oscillation frequency is 



v= -W0.12 •••6.7)- l(r 3 Hz. (2i 



Here m is the mass per area of the movable plate measured in g/cm 2 . Thus it seems 
impossible to observe even one entire oscillation, because the temperature relaxation 
takes less time. At higher temperatures the situation is slightly improved. At 245°C 
the force is 10 times and the frequency 4.2 times larger. 

4. Conclusions and outlook 

In order to discuss thermal effects, we have considered the standard setup for the 
Casimir force: two parallel, conducting square plates (side L) in a distance a. The 
plates are enclosed in a cube L x L x L with one plate at face. An ensemble of free 
photons fufilling boundary conditions at the plates resp. faces was investigated. 

In section 2 we have used Quantum Statistical Mechanics to compute the (free energy, 
entropy, internal energy) per area L 2 and the pressure. We have calculated the contri- 
butions from the space between the plates resp. outside of them separately. 
In order to keep track of ultraviolet divergences, a standard regularisation was used. 
The divergences are absent both in the entropy and in the total pressure. They survive 
both in the free and internal energy in the extensive contribution from the external 
cube (cf. Table 1). 

Separating all extensive contributions, the remainders can be written as products of a 
factor containing some power of a and a function of a dimensionless variable (8) (cf. 
(13)). These four functions are related by two linear equations, of which one (15) is 
well-known from thermodynamics. The other one (14) is connected with scaling prop- 
erties. Both formulae (Table 2) and diagrams (Fig.l) for the functions are presented. 
The results allow for accurate quantitative information, also in the (relatively narrow) 
domain between the behavior at low resp. high temperatures. 

In section 3 we have considered a situation with different temperatures between (T) 
resp. outside (T") of the plates. For T > T' the Casimir pressure is reduced by 
thermal effects and at some distance a the total pressure can vanish, so that the regions 
inside resp. outside the plates are in mechanical equilibrium. If both T and T' are 
fixed (isothermal case), this equilibrium has turned out unstable. If only T" is fixed 
and T is allowed to vary with a in such a way, that the entropy between the plates 
remains constant (adiabatic case), we have found also a region with stable equilibrium, 
if T'/T(a = 0) is small enough. The frequency of oscillations of one plate about stable 
equilibrium distances turned out too low for an observation (cf.(28)). Whether the 
(very small) restoring force (cf.(27)) can be measured, must be left open. 

It is known ||, that the Casimir force depends on the geometry of the setup. Accurate 
experimental results |], [7] have so far been obtained only for a setup consisting of a plate 
and a sphere. Temperature effects have been observed, but only the low temperature 
behavior was involved. An experiment with the original setup (as studied here), carried 
out at larger distances, would test the theory in a different domain. 
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A. Appendix 

Here we shall give a collection of formulae for the polylogarithm C, which we have used 
in computations. Some of these are also useful in Quantum Statistcal Mechanics of 
ideal Bose gasses with other constituents than photons. 

Let r, s be integers and < y < 1. The polylogarithm C r (y) can be defined by the 
series 

oo n 
Cr(v) = E ~r ( A1 ) 

and fulfills evidently the relations 

£ r (0) = 0, £ r (l) = C(r), £ r (y)>0. (A2) 

By differentiation we obtain 

y—C r (y) = A-i(y), (A3) 

dy 

f C r (y)^ = Cr+i. (A4) 

J y 

For r < 1 the function is elementary. We have 

£ 1 (y) = -ln(l-y), ( A 5) 

Co(v) = Y^- y (A6) 

Putting y = exp(— z), we obtain from (A3) 

(j-\ C r (exp(-z)) = (-l) s £ r _ s (exp(-^)). (A7) 
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From (A4) we obtain 

/ £ r (exp(— z))dz = — £ r+ i(exp(— z)). (A8) 

Combining (A7) and (A8) and using repeated partial integration, we obtain 

/s yS—n 

z s C r (exp(-z))dz = -s\ ^2 / _ y A-n»+i(exp(-z)). (A9) 



ra=0 



With r = 1 this formula can be used for the evaluation of the integrals on Gp in (3), (4). 

For partition functions we have to put z = nx and to evaluate infinite sums on n. A 
useful formula reads 



y C r (exp(-nx)) = V -, — r (A10) 

ti ^n-exp(nx)-l l ; 

The second form is obtained from the first one using the definition of the polyloga- 
rithm and summing the resulting geometrical series. It is noted, that the second form 
converges exponentially (also for negative r), as long as x > 0. For an approximate 
evaluation the series can be terminated, if x is large enough. 

The functions (10) can be obtained starting from (A10) with r — 1: 



j'( a: ) = E-3^' iV = exp(nx) - 1 > 0. (All) 

Carrying out the derivatives, we obtain 

,.. . d . . . J2, 1 / 1 + nx nx\ , . s 

k (x) = (1 - x-r-)j(x) = E-J — =- + th > 0, (A12) 



dx' v ' ^n 3 \ N N 2 



n=l 



1/1 3 2 \ 

/; \N ' N 2 ' N 3 ) 



^) = - 2 E-b + 4 + ^3 <«■ (A13) 



These forms have been used in all numerical calculations. 
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